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Abstract-This paper presents an indirect adaptive fuzzy control 
scheme for a class of unknown multi-input multi-output (MIMO) 
nonlinear dynamic systems with external disturbances. Within 
this scheme, the dynamic fuzzy logic system (DFLS) is employed 
to identify the unknown nonlinear dynamic system. The control 
law and parameter adaptation laws of the DFLS are derived 
based on Lyapunov synthesis approach. The control law is 
robustfied in H* sense to attenuate external disturbance, model 
uncertainties, and fuzzy approximation errors. It is shown that 
under appropriate assumptions, it guarantees the boundness of all 
signals in the closed-loop system and the asymptotic convergence 
to zero of tracking errors. Extensive simulation on the tracking 
control of a two-link rigid robotic manipulator verifies the 
effectiveness of the proposed algorithms. 

Keywords-MIMO nonlinear system, DFLS, Lyapunov synthesis 
approach. 

I. Introduction 

Identification and control of nonlinear systems have 
attracted a lot of attention and represents a challenging area in 
control community during the last two decades. The 
development of geometric nonlinear control theory and in 
particular feedback linearization methods, have led to great 
success in the development of controllers for nonlinear systems 
[1-2]. A key assumption in these techniques is that the 
dynamics of the nonlinear systems are exactly known. Some 
limitations of this theory appear because real systems may have 
uncertainties. Thus, to deal with uncertain nonlinear systems, 
many adaptive control approaches have been proposed. 
Adaptive control approaches are applied to systems with 
parameter uncertainties. Several results can be found in [3-6] 
and the references therein. 

The introduction of fuzzy systems led to a great success and 
provides effective approaches to handle nonlinear systems 
especially complex and ill-defined dynamic systems. Being one 
of the efficient intelligent techniques, fuzzy systems have been 
applied to the modeling and control of uncertain nonlinear 
systems. Based on the universal approximation theorem [7] 
several stable adaptive fuzzy control schemes have been 



developed for unknown single-input single-output (SISO) 
nonlinear systems [7-10], for MIMO nonlinear systems [11-13], 
and for large-scale interconnected nonlinear systems [14-16] to 
provide an effective framework for incorporating the expert 
knowledge systematically and achieve stable performance 
criterion. The stability analysis in such schemes is performed 
by using the Lyapunov synthesis approach. However, these 
adaptive fuzzy control schemes are static in nature. Motivated 
by the fact that most of physical systems are generally dynamic, 
this suggests that one may introduce some sort of dynamics to 
these static fuzzy models in order to cope with the dynamic 
nature of the physical systems. This would provide a new tool 
in the control of dynamic systems. The resulting dynamic 
structure so called the DFLS was first introduced by Lee and 
Vukovich [17-18]. They have successfully applied this concept 
to the identification of single-link robotic manipulator [17]. 
Stable identification and adaptive control based on DFLS was 
performed in [18]. This work extended to a larger class of 
SISO nonlinear systems in [19]. 

However, the previous work on DFLS is limited to only 
SISO nonlinear systems. Based on the initial results of SISO 
nonlinear systems [18-19], we intend to develop adaptive 
dynamic fuzzy control for MIMO nonlinear systems. 
Furthermore, the proposed scheme is constructed by integrating 
the feature of H°° tracking performance which can greatly 
attenuate disturbances, model uncertainties, and fuzzy 
approximation errors. 

The paper is organized as follows. A class of MIMO 
nonlinear systems and control objectives are described in 
section 2. Section 3 presents a brief description of static fuzzy 
systems and DFLS. The DFLS based-adaptive control design is 
presented in section 4. In section 5, the proposed control 
algorithm is used to control a two-link robot manipulator. 
Section 6 concludes this paper. 
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II. System Description 

In this paper, we consider a class of MIMO nonlinear 
system represented by the following set of differential 
equations: 



y 1 ( * ) = /i(*)+ZM*) I, ; +d i 



Control Objectives: Develop a feedback control law 
u(t) (based on DFLS) which ensures the boundness of both all 

variables in the closed-loop systems and the parameters of the 
DFLS, and guarantees output tracking of a specified desired 
trajectory y = [y ml , ,-, ,y m ] r • In addition for a given 
disturbance attenuation level p > , the following H°° tracking 
performance index is achieved. 



(1) 



V'WpOO+ZM*)"^ 



2 Jo 



j=i 



i_ r ■ (",-1) ■ (n„-l) n T 

where x = [y 1 ,y l , ...,y 1 , ... ,y p ,y p , ...,y p " ] 
is the overall state vector which is assumed to be available for 
measurements, u=[u 1 , ...,u ] T is the control input 

vector, y — [y 1 , . . . , y ] T is the output vector, 

D = [d 1 , . . . , d ] T denotes the external disturbance, and 

f ( x )' 9 ( x )» '» J = 1. ■ ■ ■ » P are smoom unknown nonlinear 
functions. Let us denote: 



y (n) =[yi (ni) ...y P ( " p) ] r fm = [&(*)... f p {x)f, 



T 1 1 

J e T Qe dt<-e T i (0)P i e i (0) + -h i z T z(0) 
+ U T (0)A(0) + +^p 2 f Q S T Sdt 



(3) 



G(x) 






and D 



Then, the dynamic system described by (1) can be rewritten 
in the following compact form: 



y w =F(x) + G(x)u + D 



(2) 



Throughout this paper the following assumptions are 
considered on the system (1) 

Assumption 1. The matrix G(x) is bounded away from 

singularity over compact set U CZ R n , specifically 



\\G(x)\\ =Trace(G T (x)G(x))>b 1 >0 
represents the smallest singular value of G(x). 



where 



Assumption 2. The reference trajectories y . , i =1, ••• , p is 

known bounded function of time with bounded known 
derivatives, and it is assumed to be r -times differentiable. 



where e is the error vector, 5 e L 2 [0, T\ is the combined 
disturbance and approximation error for T e |0, ooj, Q, P are 

positive matrices of proper dimensions, A is parameter 
approximation error vector, Z is identification error vector , 
and h is a design parameter. 

III. Description of DFLS 

The DFLS is composed of an ordinary fuzzy logic system 
(also referred to static fuzzy logic system) and a dynamic 
element. The basic structure of the fuzzy logic system 
considered in this paper and has been widely used in 
identification and control of nonlinear systems, is shown in 
Fig. 1, it is composed of four major components, namely, a 
fuzzification interface, a fuzzy rule base, a fuzzy inference 
engine, and a defuzzification interface. 

For MIMO fuzzy systems, the fuzzy rule base is made up of 
the following inference rule: 

R 1 : If X 1 is F/ and X 2 is F 2 ' and... X ; is F.' Then y i is 

G[ and y 2 is G' 2 and... y } is G'. (4) 

where F 1 and G 1 are fuzzy sets in R , / = 1, 2, • • • , JV ; 
i =1,2, ■■■ , n ; j = 1, 2, • • • , p . Fuzzy inference (4) can be 
decomposed and expressed as: 

R' : If X 1 is F/ and X 2 is F 2 ' and... X n is F n ' Then y } 

isGj (j=L2,-,p) 

Through center-average defuzzifier, product inference, and 
singleton fuzzifier [7], the output of the fuzzy logic system can 
be expressed as: 
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Fig. 1 MIMO Fuzzy Logic System 



^W" v N ^TT" T^T (5) 



where y is the centre of the fuzzy set G at which 
jU G achieves its maximum value, and we assume that 



Using the definition in (8), equation (9) can be written in the 
following compact form: 



y = -ay + Y T 0(x) 



(10) 



/4(y ! )=i 



Eq. (5) can be written as: 

y j (x) = Y j T t(x), j=i,2,--,p 



(6) 



where y = [y 1 ,---, 9 ] is me output of MIMO DFLS, 
a = diag[a l ,■•■, a ] is a positive constant matrix. 

The DFLS described by (10) was shown to possess universal 
approximating capabilities to a large class of nonlinear 
dynamic systems [17]. 

Let Z e R p a vector defined as: 



where Y —[y 1 ,---, y N ] T is a vector of adjustable 
parameters, and ^(x) = [^ , • • •, (/>„ ] is a regression vector 

with each (pi variable defined as a fuzzy basis function (FBF) 
[7] as: 



y 



("i-i) 
i 



y { r ] 



(11) 



The system (1) can be written as: 



^=^7 



nx.w 



H(ir^w) 



(7) 



Define Y = [Y l , • ■ • , Y ] /x as a matrix of adjustable 
parameters. MIMO fuzzy system can be expressed as: 









(12) 



y=Y m 



j-g-v Or equivalently, equation (12) can be written in the compact 
form as: 



The MIMO DFLS shown in Fig. 2 can now be described by the 
following differential equations: 



z = F(x) + G(x)u+D 



(13) 



9i =- a i9i +Y i T <f>(x) 

y P =-a P y P +Y p T 0M 



According to the universal approximation theorem [17], the 
following DFLS can be used to identify the unknown MIMO 
(9) nonlinear system (12). 



y = -ay + Y T <fi(x,u) 



(14) 
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Fig. 2 MIMO Dynamic Fuzzy Logic System 

Our objective now is to develop an appropriate control law for For reference trajectories to be asymptotically tracked we 

input u in (1), and an adaptation law for the parameter matrix choose: 

Y of the DFLS, such that the closed loop system is stable in 

the sense that the tracking errors e , = y ■ — y , , j — 1, • ■ • , p , 

are uniformly bounded and at the same time, the identification 

errors and identifier parameters are also uniformly bounded. 

IV. DFLS Based Adaptive Control 

In this section, we develop an adaptive control scheme for 
system (1) based on DFLS based identification. 

Consider the system (1), for the given reference 

trajectories y m =[y ml , •■• , y mp ] ■ Let us define the tracking 

errors as: 



v = v (n ° +k e (ni ' + --- + k e 

v l J ml T * v lr 1 c l T,v ll c l 

(n p ) 



(18) 



+ k P r p e P 



V = V 

p J mp 

Substitute (16) in (1) yields 



( " p_1) +-+k e 
T pi p 



,("i) 



+ /c l „e 1 ( " l " l) + --- H-k^e, =0 



v lr, "1 



11"! 



(19) 



e i y m \ y^ 

p J mp J p 



( n n) i ( n „-l) 7 

e p ^V P e p ^^p^p 







(15) 



If the coefficients /C- are chosen such that all polynomials in 

(19) are Hurwitz stable, then we can conclude that 
lim e ; (t) = which is a main objective of control. 



Denote e = [^ , • • • ,e p ] T then e = y m - y 

If the system (1) is known i.e. F{x) , G(x) are known and Consider the MIMO DFLS in the form of (9): 

z 1 =-a 1 z 1 +Y 1 T 0(x,u) 



D = , the feedback law: 

fln(x) ■•• 9i P ( x ) 

9 P M)--- 9pp(x) 
Yields the linearized systems 



-1 


f 


"AM" 


+ 


1^ 

v l 




V 


.f P W. 




- V JJ 



(20) 



(16) 



2 p =-« p z p +Y p VO,u) 






(«i) 



which can be used to identify the unknown MIMO nonlinear 
system (12). Define identification errors as 
z=z-z., i = l,2,-p- Our objective is to develop an 

appropriate control law for input u t and an adaptive law for the 

(17) identifier parameters Y- such that closed-loop system is stable. 
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The expression for Z f . in equation (12) can be written as 

z 1 = -a 1 z 1 +Y 1 T <fr(x, u)-r 1 (x,u,^,7 1 ) 
2 P =-a p z p +Yp0(x,u)-r p (x,u,(f),Y p ) 



In this situation we propose the following control law which is 
based on DFLS: 



(21) 



where r,(x, u,^,^) represents the static modeling error of the 
DFLS identifier and can be expressed as: 



9n(x) ••• 9ip(*j 
SpiOO'" 9 PP ( X ) 



a l 





o" 


z l 


















«p_ 


_v 



Y p T <f>(x,0) 


+ 


_V 


+ 





(27) 



r i (x,u,^,Y / ) = -a 1 z 1 +Y / r ^(x,u) 

-f,00-£Mx)«,-d, 



(22) 



Or equivalently, equation (27) can be written in the compact 
form: 

u = G(x)~ 1 [az-Y T 0(x,O) + v + u r ] (28) 



By Lemma 1 in [17], there is exist an optimal parameter where G(x) is a static fuzzy logic estimation of G(x) , 
vectors 

^(x,0) = <j>{x, u) L =0 , and U r is a robust compensator which 



y.*=min{y. : Y. <M.} 



which minimize the static modeling error, V i , such that 



(23) 



defined as 



1 T 

u ri =—B?P : e : 

x. 



I I —I 



(29) 



sup 

(x,u)eQ 



r t (x,u,^X) 



<m: 



(24) where X i , P t are the solution of the following Riccati-like 
equation. 



where M- and M r are positive design constants. In the 

following, we develop an adaptive law for Y i . Replace Y t by 
Y t in (21) results in 



AfP i+ PA -Q, 



r 2 ^ 



z r =-a l z 1 +Y 1 * ^(x,u)-r 1 (x,u,^,Y r 1 *) 
i P =-^ p z p +Y* T ^){x,u)-r p (x,u4,Y*) 

Subtracting (25) from (20) yields 

Zj --a 1 z 1 +/S T l (j){x,u) + r l {x,u,(j),Y*) 
Zp =~ a pZp + AT p 0( x ' u ) + r P ( < x ' u '0'Yp) 

where A. —Y ( — Y t * is the parameter estimation error. 



K p 

Using (27), we can rewrite (12) as 



PB ; BjP=0 (30) 



( ( ( I 



(25) 



f P (x) 



(g(x) + G(x)-G(x))g- 1 (x) 



(26) 



a x 
'•. 
« 



Y^ixfi) 



Ylftxfl) 



|— — | 




— — 1 


v l 


+ 


"n 


_v 




_u r2 _ 



(31) 
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Using (17), (18), (31), and after straight forward manipulations, 
it can be easily to obtain 



The control gain functions g .. (x) can be approximated by a 
static fuzzy logic system (6). It follows that 



,("i-i) 



e w + k e ( " 1 ^ + ---+k e 



( n i>) 7 ("n _ l) 7 

e p +k P r/p P +---+^ P ie p 



gM = gJx) + K ii =YU(x) + K il 



(36) 



f 1 (x) + a 1 z 1 -Y 1 V(x,0) 



f p (x) + a p z p -Y p T </>(x,0) 



where Y is an adjustable parameter vector, <fi(x) is a fuzzy 
basis function and K^ is the fuzzy approximation error. 



"l 


+ 


""rl" 


+ 


"di" 


_v 




_ U r2_ 




_ d pJ 



+ (Ct(XJ Ct(XJ) (-32) Define an optimal parameter estimates Y.- such that it 

minimizes the approximation error. Therefore, we can write 

9yW = yy*VW + 4 (37) 



It is clear from (32) that U ri can attenuate the external 
disturbance and fuzzy approximation errors. 

Equation (32) can be written as: 

e^+k ir e^ + --+k iA = f i (x) + a i z i -Y i T ^(x,0) + 
£(M*)-Mx)K-u ri +d,, '=1>2,--P (33) 

The state-space form for (33) can be written as: 

7,(x) + a 1 z,-^(x,0)' 



where /f ; - is the minimum approximation error. Using (36) 
and (37), equation (35) can be written as 



e t =A i e i +B i u n +B i 



A^(x,0)-£A^(x)u. 



(38) 



+ B,w, 



e, = Ae, + B,-« ri + B i 



t{g l] (x)-g i] (x)} i] 



;=i 



where Ay = F y T - Y y * and w, = *J + r, (x, u,0, Y,* ) . 
The adaptive laws are chosen as 

Y t =-jj i h i z i <f>(x > u) + jj i <t>(x > 0)B?P i e i 



(39) 



where 



v i(n,-l) i(n,-2) 



B, 



Bd. 

I i 

(34) 



and 



Y.. =-T lij e T i P i B i <f>(x)u j 



(40) 



Using (22) with u = , we can write 
e, = Ae, + B,u . + B ; 



i 11 in 1 



- A^(x,0) + £( ffj/ (x) - g y (x)) Uj . - r ; (x,u,0, 



(35) 



Theorem 1. Consider an unknown MIMO nonlinear dynamic 
system (1) which is controlled by (27) and to be identified by 
the DFLS (20) by adjusting the parameter vectors Y T and 

Y. T with the adaptive laws (39) and (40) respectively, then the 
closed loop system possess the following properties 

(i) All signals in the closed-loop system are uniformly 
bounded. 

(ii) For a given disturbance attenuation level, the proposed 
tracking performance index (3) is achieved. 

Proof. Choose a Lyapunov function as 



V=V, 



V 

y p> 



(41) 



JCET Vol.1 No.l July 2011 PP.23-33 www.ijcet.org ©World Academic Publishing 



28 



Journal of Control Engineering and Technology (JCET) 



1 1 1 p 1 

2"' '"' 2 2/7,. ^2/7, 



A^A, (42) 



Differentiating V, and Vj 



v = y 1+ ...+v 



(43) 



Substituting (26) and (38) in (43) after differentiation, we can 
obtain equation (44) at the bottom of this page. 

Using the adaptive laws (39) and (40), equation (44) can be 
simplified into 



2 Pi 



2~2 ^2 



— — — H r ' 

2p 2 2 



(47) 



+ ^(w l BfP i e i +e T i P i B i w i ) 



Note that the third and fourth terms are negative in (47). 
1 



Vj < -e*Qe t -e T : P : B : Bje ; 



2p 



2—i i i i —i 



( 



V> <e 



A T P,+PA -je*P t B t Bfe t 



e : 



(45) 



"*--!P 



V 



z t 2 +?-r t \x,u,tX) ( 48 ) 



J 



■a,h,z, ! -h,z, r,(x,u,^,Yi-)-i(w,B, r P,e, + efi>B,w,) 






+ ^{w i BfP i e i +e]P i B i w i ) 



Using the following triangular inequality for the third term in 
(45) 



The third term in (48) can be made negative by choosing 
h,. <2a iP 2 . 



h?r P 2 



h,z t r^u,^ )£-i-L- + l_r/(x,u,£r, ) (46) 

2p 2 



v^-£Q^-\ 



1 r 

V/ 7 



P-( ri 2 (x,u,<pX) + ™?) 



(49) 



: T n _ , 1 j n • , ,. ~ ~ , 1 ;r a , V ^ iT 



Vi =^-e;P i e 1 . +^e;P,e ; . +h i z i z t +— A^A,. +T — A^A f 
2 2 77,. ^/7, 7 



1 p 

ef AfP,.e,. -— efP^Bfe,. -^x,0) T A i BfP t e i -^^x)A ij BTe i u J +w i BfP i e i 



p 



1 
+ efP ; -A- e t -— efP I B,B, T e i -efP,B ! A^(x,0)-^efP I .B,.A^(x)u j +efP f B,w,. 



H-hjZ,. 



-a i z l l +A i <f>(x,u) + r i (x,u,<f>,Y i ) 



1 - T 
7,- 



' + ^7 



1 ,-^r 



V A (/ 



= e r A T P+PA e T PBB T e 



A 



e z . -a,h,.z, 2 +— (y, t +77,.h ; .z^ T (x,u)-77^(x,0)B, T P,e J . k 
77,. 



+ 



-7 y 



£— (y,/ +/7 ij ef P ; B^ T (x)u ; .)\ l7 +h,z i r i (x,u,M-V^(w I -BJ r i>.e I . H-ef^-B^ 



(44) 
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Since 



1 T 
— e i P i B-pw i 

P 



, from (49) we obtain 



Let 

Q = diag[Q 1 ,-,Q P ] , P = diag[P 1 ,---,P p ] 

It t y ~ r~ ~ ir 
?i.-»?pJ • z = L z i'--' z P J > 



^^-e-Qe,+^ 2 ^ 2 



(50) 



J_A r ... — A T — A T ••• — A T 
7i 7 P 7n 7 PP 



and 



where J. 2 = r ; . 2 (x,u,^, Yj*) + w 2 . 

After some straightforward manipulations, we can deduce 

V.Z-cK+Mt (51) 



s=W-,s p ] 



where 



c, = min</l. 



1 1 



7, 7, 5 



with A = inf ^ in (Q,0 and 
SU bA max (Q,) 



Then from (54) we obtain 

jVQe tft<-e r (0)Pe(0) + -hz r z(0) 

+ ^A r ^)A(0) + + VjV£dt 



1 r r 
2 



(55) 



V < -cV + ju 



(52) 



1 p 



From (51) and (43) 



where 






This implies that all signals in the closed loop system is 
bounded. Thus the control objective (i) is realized. 

Integrating (50) from t = to t = T , we have 



Thus the control objective (ii) is achieved and the proof of 
Theorem 1 is completed. 

V. Simulation Results 

To demonstrate the effectiveness of the proposed scheme, a 
simulation is performed for the tracking control of a two-link 
rigid robot manipulator moving in horizontal plane. The 
dynamics of the robot are described by the following 
differential equation: 



(56) 



M u M 12 ~ 


~<Zi~ 










M 21 M 22 


32 _ 






- Mi ~ hq 2 ~ Mi 


Qi 




?1 


-H 







<J2 




T 2 



where 



[ T e T Qedt<V(0)-V(T) + -p 2 [ T S 2 dt (53) Af u = a 1 + 2a 3 cos(q 2 ) + 2a 4 sin(q 2 ) 

Jo -' ' ' ' v 9 Jo ' 



1 f T 

2 



Since V (T) > we can write (53) as follows: 



J o r ^Qe,dt<V I .(0) + i / 7 2 [ r j 2 d t 



1 r T 
2 



^■(0)P i e,(0) + ^ I z j 2 (0) 



1 



(54) 



+ -^A^(0)A,(0) 
27, 






M 22 =o 2 

M 21 = M 12 = a 2 +a 3 cos(q 2 ) + a 4 sin(q 2 ) 

ft = a 3 sin(q 2 ) - a 4 cos(q 2 ) 

with 

a, =L+ m,/ 2 , +L+ ml* + ml? 

1 1 lcl e e ce el 

a 3 = m e l 1 l a cosS e 
a 4 = m e l 1 l ce sinS e 

In the simulation, the following parameter values are used: 



30 
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^=0.1 / cl =0.5 1^=1.0 ^=0.12 
I ce = 0.6 S e = 0.6 m e = 2.0 I e = 0.25 

Since the inertia matrix M is positive definite, the system can 
be written as: 



Ri 
q. 



+ 



M u 


M 12 " 


-i 


-hq 2 


M 21 


M 22 _ 




r hq t 


~M U 


M 12 " 


-i 


Ti 




[_M 21 


M 22 




J 2. 





hq 2 -hq l 




<j 2 



(57) 



// p5 = exp(-5(x,-0.5) 2 ) , 

y" F(6 =exp(-5(x,. -I) 2 ) , 

/^ = 1/(1- exp(-5(x ( -1.5))). 

For given 

Q, = diag[10X0] , p = 0.1,0.2, A = 0.01,0.02 . 

Solving Riccati equation yields 



P. 



15 5 

5 5 



Let x 1 = q 1 , 


*2 =C ll> 


x 3 = q 2 , 


y 2 = x 3 


G( 


F(x) = M 1 


~-hq 2 


-hq 2 - h Ri 




G(x) = M~ 1 
Qi 

Then, the dynamics of the two-link robotic manipulator can be 
expressed as: 

y = F(x) + G(x)u (58) 

The control objective is to force the system outputs q 1 and q 2 
to track the desired trajectories: 

y ml =0.2sin(t), y m2 = 0.2sin(t) 

For the DFLS, seven membership functions are selected as 
follows 

ft n = 1/(1+ exp(-5(x,.+ 1.5))), 
M F , = exp(-5(x,+l) 2 ), 
// F3 =exp(-5(x.+0.5) 2 ), 

ju , = exp(-5x, 2 ) 



Choosing 

a 1 =9,a 2 = 5, \ = 1000, h 2 = 200, 

k u = ^21 = 4 ' Kl = ^22 = 10 > ^i =Vi =1» 
7 n =7 22 =0.1, 7 U =0.11, 7 21 = 0.21 
Choose initial condition as 

Xj(0) = x 3 (0) = 0.1, x 2 (0) = x 4 (0) = , 

and the initial conditions for the adaptive parameters are 
chosen to be zero. Figs. 3-6 show the simulation results. 

VI. Conclusions 

In this paper an adaptive fuzzy control scheme is 
developed based on DFLS for a class of uncertain nonlinear 
MIMO systems. DFLS is used to identify the unknown 
nonlinear system as a whole and an adaptive fuzzy controller is 
developed. The fuzzy control law is robustfied by an H 00 
compensator to attenuate disturbances and fuzzy approximation 
errors. The proposed approach guarantees that all signals in the 
closed loop system are uniformly bounded. Simulation results 
show the effectiveness of the proposed scheme. 
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Fig. 3. Tracking curves of qi 
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Fig. 4. Tracking curves of q 2 
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Fig. 5. The control input Ui 
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Fig. 6. The control input u 2 
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